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Abstract We present a new method for the analysis of Abell 1835 observed by XMM-Newton. The method 
is a combination of the Direct Demodulation technique and deprojection. We eliminate the effects of the 
point spread function (PSF) with the Direct Demodulation technique. We then use a traditional deprojection 
technique to study the properties of Abell 1835. Compared to that of deprojection method only, the central 
electron density derived from this method increases by 30 percent, while the temperature profile is similar. 
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1 Introduction 

Being the largest and most massive celestial bod- 
ies in the universe, galaxy clusters are the ideal labo- 
ratories for studying large-scale structure of the Uni- 
verse as well as dark matter, through which the cos- 
mological model^ could be tested. The center of 
galaxy cluster is extremely complex, and the observa- 
tion is susceptible to interference. The XMM-Newton 
observatory is an ideal instrument for observations of 
galaxy clusters for its excellent spatial and spectral 
resolution. EPIC pn and MOS are the main instru- 
ments onboard XMM-Newton designed for X-ray de- 
tection. However the PSF of pn and MOS is not ig- 
norable, which makes the spatial resolution worse and 
causes the analysis at the central cluster inaccurate. 
XSPEC could eliminate the effects of PSF through 
fit-ting several spectra parameters . However, this 
method can not be used with deprojection simul- 
taneously, which is an open problem for many years 
in the analysis of clusters, and the fit models here do 
not consider the energy dependence of the PSF. 

The direct demodulation (DD) is an effective 
method for image restoring . Through Lucy 
iterations 1-1 with constraints of background - , the re- 
stored image has better spatial resolution, helping us 
learning of galaxy clusters greatly. The DD combined 
with deprojection is a new method which can do PSF 
correction and deprojection simultaneously without 
empirical model completely. Comparing to the above 
method, the DD is more effective for the sensitivity 



and spatial resolution, even if the statistic and signal- 
to-noise ratio is bad. 

Abell 1835 (z=0.2523) has a relaxed structure in 
which a large cooling flow has been found by ROSAT, 
ASCA and XMM-Newton M _ xhc i um i n0 sity 
of Abell 1835 is sufficient to DD, and its relaxed and 
symmetry structure is suitable for deprojection. In 
this paper we restore pn images in different energy 
channel of Abell 1835 (observation ID 0098010101); 
and then we employ a traditional deprojection tech- 
nique applied in several clusters successfully . Fi- 
nally, the spectra with-out PSF and projection ef- 
fects are obtained. This is the first time to deproject 
a galaxy cluster considering PSF. 

The structure of this paper is as follows: section 
2 describes PSF, DD and processes the data of Abell 
1835 with DD; section 3 presents the deprojection 
method, the fit results such as the temperature pro- 
file and the electron density, and the error estimation; 
then we conclude the results in section 4. 

2 Direct Demodulation 

We can use the Lucy iteration'-' (eq. 1) with back- 
ground constraints to restore the sub-images of Abell 
1835 in 34 different energy bands. 

i' 

(1) 

where /' r+1 '(i) can be derived from the initial 
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values,/' - 1 (i), with several iterations and background 
constraint. 

The average value of each energy band is taken for 
the calculation of PSF, which is a function of energy 
and angle. The steps to process a sub-image are as 
follows: (a) Select a region of 100x100 with the cen- 
ter (263, 279) as the DD region, about 500500 arcsec 
square (1 bin = 5arcsec); (b) Find the positions of 
the gaps between pn CCDs; (c) Select a circle, whose 
center is the most luminous in the sub-image and 
the radius is 120 pixels; calculate rest region with- 
out CCD gaps as back-ground; (d) Iteration with the 
above background; (e) Combine 90x90 of DD region 
with the rest of original sub-image into an intact sub- 
image, which can avoid the edge error of DD region. 

Because the DD processes the iteration with phys- 
ical limit, its error is hard to obtain through the di- 
rect error estimation. The error depends on the PSF, 
background, shape of source and iteration times. We 
find that the 10 times iteration is the best choice 
for DD. Obvious dispersion will appear in the im- 
age when the iteration time is more than 10, while 
the result is insufficiently stable when less than 10. 
For estimating the errors of DD, we create 20 images 
with the Poisson distribution by convolving the result 
/ and the PSF of each energy band; DD the 20 im- 
ages one by one and get 20 different values of / which 
is f'jU = 1, •••,20). We use the standard deviations of 
/' as the errors of /. 

We take the 0.3-0.4keV sub-image as an example. 
The original and the PSF-corrected images are shown 
in Fig.l. The counts in each ring are calculated be- 
fore and after DD. The restored image shows a much 
higher counts in the central region (Fig. 2). 
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Table 1. The fit results of the spectra. Tl and 
T2 are the temperatures in the double temper- 
ature model. Al and A2 are the abundances 

in the double temperature model. The error Fi S- 1 - The to P is the original image; the bot- 

bars are at the 90% confident level.. tom is the ima S e after PSF correction. 
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Fig. 2. The average counts in each ring before 
and after DD. Triangles and squares represent 
the original counts and the counts after PSF 
correction, respectively. 

3 Deprojection and spectral analysis 

3.1 Deprojection 

Assuming the structure of Abell 1835 is spheri- 
cally symmetric and the spectra of each ring for the 
deprojection are the same, then the deprojected spec- 
tra could be calculated through eliminating the con- 
tribution from the outer rings for all spectrum compo- 
nents. We divide the sub-image into 8 annular areas, 
whose center is located in the most luminous pixel, 
and use the utmost ring (6' — 8.33') to obtain the X- 
ray background. However, the signal in the seventh 
ring is almost background, so we only consider the 
inner sixth rings (^6'). 

3.2 Spectra analysis 

Since different version of SAS and RMF of pn 
could derive different results, so we use SAS 8.0 and 
XSPEC 12.4.0 in this paper all the time. We ana- 
lyze the DD and deprojected spectra of pn data in 
XSPEC with a single and a double temperature ther- 
mal plasma model (eq.2 and 3) as follows: 

Model 1 = Wabs(n H ) x Mekal(T,z,A,norm) (2) 

where Wabs is a photoelectric absorption model^ 
and Mekal is a single temperature plasma emission 
model 1111 . 

Model 2 — Wabs(n H ) x (Mekal(T 1 ,z,A 1 ,norm 1 ) 
+Mekal(T 2 , z, A 2 , norm?)) 

' (3) 

The latter is the double temperature model implying 
there are two components in the central region with 
different temperatures. 



Because the outer region (3.33' — 6.0') contains 
pixels out of DD region, we ignore this ring and its 
adjacent ring (2.25' — 3.33'). We fit the three rings 
(0' — 1.5') of the central region with a double tem- 
perature model with free abundance. It was found 
that the region (0.75' — 1.5') can be also fitted well 
by a single temperature model like the outer region 
(1.5'-2.25'). 

As shown in Table 1, the normalization of the high 
temperature component in the innermost region is 
about 60 percent higher than that of deprojection 
only. It means that the central electron density in- 
creases by about 30 percent, implying a much larger 
gradient of the gas density in the central region of 
Abell 1835. However the resultant temperature pro- 
file is more or less the same as that of the deprojec- 
tion only, except for the innermost region where the 
temperature of the high component increases signifi- 
cantly. 

Since the profile of PSF is dependent on the pho- 
ton energy (Fig. 3), the temperature of the restored 
spectra could be different from the original. The 
PSF of high energy is steeper than that of low energy 
around the several adjacent bins, especially in the 
inner 2 bins. In addition, the changes of the normal- 
ization of the restored spectra will affect the spectra 
of its inner rings during the deprojection process. 
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Fig. 3. The PSF of pn. The heavy solid line, 
the dotted line and the light solid line repre- 
sent the PSF of pn of 0.3keV, 3.8keV, 7.8keV, 
respectively (1 bin = 5 arcsec). 
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Fig. 4. Spectra fitted by the double tempera- 
ture model with free abundances. The top is 
for the ring of 0' — 0.25'; the bottom is for the 
ring of 0.25' -0.75'. 

The spectra of the two innermost rings (0' — 
0.25', 0.25'— 0.75') fitted by the two temperature mod- 
els are shown in Fig. 4. 

4 Conclusion 

We have applied a combined DD and deprojec- 
tion technique to analyze the pn data of Abell 1835 
observed by XMM-Newton observatory. The results 
show that it is reasonable to analyze the cluster data 
with the combined DD and deprojection technique. 
Comparing to the analysis by deprojection only, we 
derive a similar temperature profile but a higher cen- 
tral electron density. This indicates that the effects of 
PSF may be not important for the temperature profile 
but it should be substantial for the determination of 
the central gas density as well as some other physical 
properties in the central region, e.g., the cooling-flow 
rate. 

The authors are grateful to the anonymous refer- 
ees for their insightful suggestions. 
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